High-spin yrast structure of "^^S suggested by 
symmetry-unrestricted, cranked Hartree-Fock 
-Q^- ■ calculations 

Masayuki Yamagami and Kenichi Matsuyanagi 
^! Department of Physics, Graduate School of Science, 

^ ■ Kyoto University, Kitashirakawa, Kyoto 606-8502, Japan 

^ ; 19 August 1999 

> 
O 

O ' Abstract 

00 ; 

\ High-spin yrast structure of ^^S is investigated by means of the 

2^ ' cranked Skyrme-Hartree-Fock method in the three-dimensional Cartesian- 

mesh representation without imposing restrictions on spatial symme- 
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the hyperdeformed-like configurations takes place on the yrast line at 
angular momentum / ~ 24, which corresponds to the "band termina- 
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yrast states in the range 5 < / < 13. 
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1 Introduction 



Since the discovery of the superdeformed(SD) rotational band in ^^^Dy, about 
two hundreds SD bands have been found in various mass (A=60, 80, 130, 
150, 190) regions |jl|-§[. It turned out that every regions of superdeformation 
have their own characteristics and offer a number of interesting questions; 
investigations of them have been significantly enlarging and deepening our 
understanding of nuclear structure. Yet, the doubly magic SD band in ^^S, 
which has been expected quite a long time [|7 HlO|] remains unexplored, and will 
become a great challenge in the coming years 0. Exploration of the SD band 
in ^^S will certainly give a strong impact toward understanding the possible 
connection between the SD structure and the molecular resonance phenomena 
associated with the ^^O +^^0 configurations (see, e.g. |]lT],|T2| for reviews). 
More generally speaking, the nucleus ^^S seems to be situated in a key position 
in the investigation of possible relationships (such as discussed in [ITsHr^) 
between the SD states systematically observed in heavy nuclei and the cluster 
structures widely observed in light nuclei (see, e.g. [l^ for a review). Thus, 
excited states in ^^S have been theoretically studied by Nilsson-Strutinsky 
approaches |]7|- |T0|] , self consistent mean-field approaches |17,18|, spherical shell- 
models |T^,^, and cluster-structure and molecular- resonance points of view 



2TH25|. 

The aim of this paper is to study the high-spin yrast structure of ^^S 



from 



the point of view of exploring exotic shapes in nuclear high-spin states by means 
of the cranked Hartree-Fock (HF) method with the use of the Skyrme forces |26, 
27[| , which is called "the cranked SHF method". One of the recent advances 
in nuclear structure theory is that it has become possible to carry out the 
HF calculation in the three-dimensional (3D) Cartesian-mesh representation 
^,|32|. This approach has been extended 0,|3^,|3^ to a rotating frame 



by introducing the cranking term and applied to the high-spin yrast states of 
^^S in Ref. [1^ with the use of the BKN interaction ||31[. In these cranked 
HF calculations, however, parity and signature symmetries are assumed for 
the intrinsic wave functions in order to simplify the calculation. We refer 
an excellent review by Aberg, Flocard and Nazarewicz P] for an overview 
of studies on nuclear shapes in terms of various kinds of mean-field theory, 
especially other than the cranked SHF approach. 

Recently, we constructed a new computer code for the cranked SHF cal- 
culation based on the 3D Cartesian-mesh representation, which provides a 
powerful tool for exploring exotic shapes (breaking both axial and reflection 
symmetries in the intrinsic states) at high spin in unstable nuclei as well as in 
stable nuclei. As a first application of this new code, we investigated the high- 
spin yrast structure of ^^S, and found |^ that 1) a drastic structure change 
may occur above angular momentum / ~ 24 in the yrast line, and 2) non- 
axial octupole deformations of the Y31 type arise in the yrast line in the range 
5 < / < 13. The present paper is intended to give a more detailed account of 
this work. Quite recently, Molique, Dobaczewski and Dudek pS] investigated 
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several SD configurations in ^^S (not restricted to the yrast states) as well 
as in neighboring odd-A nuclei by means of the cranked SHF method with 
the SLy4 force |37] in the harmonic oscillator basis. On the other hand, they 
did not discuss the yrast states above J ~ 24 as well as non-axial octupole 
deformations, which are the major subjects of this paper. 

After a brief account of the cranked SHF calculational method in Section 2, 
an overview of the obtained yrast line for '^^S is given in Section 3. In Section 4, 
we discuss properties of the high-spin limit of the SD band, paying special at- 
tention to a band-crossing phenomenon associated with the level crossing with 
the rotation-aligned [440] | level. The result of the cranked SHF calculation 
is compared in Section 5 with that of the cranked harmonic oscillator (CHO) 
model calculation. In Section 6, effects of the rotation-induced, time-odd com- 
ponents in the selfconsistent mean field on the properties of the SD band are 
briefly discussed. In Section 7, we discuss about the Y^i deformed solutions 
of the cranked SHF equations, which constitute the yrast line in the range 
5 < / < 13. Although, at the present time, experimental data directly com- 
parable with our theoretical calculations seem to be unavailable, we briefly 
remark in Section 8 on some recent experimental references. Conclusion is 
given in Section 9. 



2 Cranked SHF calculation 

The cranked HF equation for a system uniformly rotating about the x-axis is 
given by 



6 <H 



IjJr. 



J,; >= 0, 



(1) 



where Wrot and Jx mean the rotational frequency and the x-component of 
angular momentum, and the bracket denotes the expectation value with re- 
spect to a Slater determinantal state. We solve the cranked HF equation for 
a Hamiltonian of the Skyrme type by means of the imaginary-time evolution 
technique pSf in the 3D Cartesian-mesh representation. We adopt the stan- 
dard algorithm p5|-pil|,[5^ in the numerical calculation, but completely remove 
various restrictions on spatial symmetries. Namely, we basically use the pro- 
cedure developed and applied to the yrast line of ^"^Mg by Bonche, Flocard 

m 



and Heenen 



except that the parity and the signature symmetries are not 
imposed on the individual wave functions. In this connection, we mention that 
a similar HF code (with parity projection but without the cranking term) was 
constructed by Takami et al. and successfully applied to the description 
of cluster structures in light nuclei, ^Be, ^^C, ^^O and ^''Ne. The same code 
(but without parity projection) was recently used to explore exotic shapes in 
proton-rich N Z nuclei in the ^°Zr region [^,^, and tetrahedral and tri- 
angular shapes are suggested to appear near the ground states of some nuclei 
in this region. In Refs. [| 



9|,^, the pairing correlations were taken into 



account in the BCS approximation. In the present calculation, we neglect the 
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pairing, since they are not expected to play an important role at high-spin 
states in ^^S. 

When we allow for the simultaneous breaking of both reflection and axial 
symmetries, it is crucial to accurately fulfill the center-of-mass condition 

AAA 

<Y,Xi>=<Y,yi>=<Y,Zi>=Q, (2) 

i=l i=l 1=1 

and the principal-axis condition 

AAA 

< J2 ^iVi > = < Vi^i > = < ^i^i >= i^) 
i=l 1=1 i=l 

For this purpose we examined several techniques and confirmed that the 



constrained HF procedure with quadrupole constraints [42| works well. Thus, 



we replace the "Routhian" R = H — u^otJx in eq. (1) with 

3 A 3 A 

R' = R-Yf^k < Y(^k)i >^ - XI f^k,k' < >^ • (4) 

fc=l i=l k<k' i=l 

In numerical calculations, we confirmed that the constraints (2) and (3) are 
fulfilled to the order O(10~^^) with values of the parameters ~ O(IO^) and 
~ 0(1). We solved these equations inside the sphere with radius R=8 [fm] 
and mesh size h=l [fm], starting with various initial configurations. The 11- 
point formula was used as the difference formula for the Laplacian operator. 
As usual, the angular momentum is evaluated as Ih =< >■ 

In this paper, we use the standard SlII and SkM* forces. With the use of the 
SIII force Tajima et al. |Q carried out a systematic SHF+BCS calculation 



for the ground-state quadrupole deformations of nuclei in a wide area of nuclear 
chart. They have carefully examined the possible error due to the use of the 
mesh size /i=l[fm] and found that the deformation energies obtained with this 
mesh size are quite accurate. On the other hand, the SkM* force p7[ 



was 



designed to accurately describe properties at large deformations like fission 
barriers, so that it may be suited for the description of super deformations |^ . 



In recent years, several newer versions of the Skyrme forces have been proposed 
(see, e.g. Ref. [^) in order to improve isospin properties of the Skyrme 
forces. Although the major purpose of them is to better describing neutron- 
rich unstable nuclei, it will also be interesting to employ such versions to 
examine the dependence of the results reported in this paper on the effective 
interactions adopted. We defer such a more extensive calculation for the future. 



3 Structure of the yrast line 

The calculated yrast line is displayed in Fig. 1, and angular momenta and 
deformations of the yrast states are drawn as functions of rotational frequency 
in Figs. 2 and 3. In these and succeeding figures, the calculation were done in 
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step of Acjrot = 0.2 MeV/^, and the calculated points (indicated by symbols) 
are smoothly interpolated by lines. The quadrupole deformation parameters 
P2 and 7 are defined as 

Ayr ^ A 

/?2 COS7 = ^ < E r-Y2o{ei, ct>,)> / <Y.tI>, (5) 

^ i=i i=i 

/32sin7 = -- < Y."^]{Y^2{Q^A.')^ym.^■^)> I <Y.^\> 

It is seen from Figs. 1-3 that the results of calculation with the SIII and 
SkM* forces are quite similar: For both cases, the expected SD band becomes 
the yrast for J > 14, and it exhibits a singular behavior at about / ~ 24. 
As we shall discuss in the next section, this is due to a band crossing associ- 
ated with the rotation-aligned [440] | level, and we call the yrast states above 
/ ~ 24 "hyperdeformed (HD)-like configuration" in order to distinguish them 
from the SD configuration. This configuration becomes unstable against fis- 
sion for / > 34. In addition to the SD and HD-like configurations mentioned 
above, we found that the yrast states with 5 < / < 13 possess an appreciable 
amount of non-axial octupole deformation of the Y31 type, so that we call, for 
convenience, this region of the yrast line "Isi band", although, as discussed in 
Section 7, some cautions are necessary in using this terminology. 

Thus, the calculated yrast line can be roughly divided into the following 
four regions: 

1) / < 4, weakly prolate region, 

2) 5 < / < 13, Isi deformed region, 

3) 14 < / < 24, SD region, 

4) 26 < / < 32, HD-like region. 

Below we first discuss the properties of the high-spin limit of the SD band, 
and later about the I31 band. The lowest-spin region will be touched upon in 
Section 7 briefly. 



4 High-spin limit of the SD band 

As we saw in Figs. 1-3, the solutions of the cranked SHF equations correspond- 
ing to the yrast SD configuration are obtained from J = to about / = 22. 
Figure 4 shows the potential energy function for the SD state at / = 0, 



evaluated by means of the constrained HF procedure with the quadratic 
constraint on the mass-quadrupole moment. We see that the excitation energy 
of the SD state at / = is about 12 MeV. 

A particularly interesting point is the behavior of the SD band in the 
high-spin limit: It is clearly seen in Figs. 2 and 3 that a jump occurs both 
in the angular momentum / and the quadrupole deformation /32 at cUrot — 2.9 
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MeV/^. At this point, / jumps from about 22 to 26, and P2 suddenly increases 
from about 0.6 to 0.7. Such a discontinuity is well known |Q to occur in the 
description of the band crossing phenomena within a standard framework of 
the cranked mean-field approach. The point is more clearly seen in Fig. 5 as a 
singular behavior of the dynamical moment of inertia J'^'^^ = dl /dumt near the 
band crossing point. (Other properties of J/^^-* will be discussed in the next 
section.) 

Figure 6 displays the shape evolution of the SD band as a function of 
angular momentum in the (/32,7) plane: With increasing angular momentum, 
small triaxial deformations gradually set in and at / ~ 24 the shape exhibits 
a striking "back-bending" toward larger prolate deformations. Evidently, this 
is due to the band crossing mentioned above. Such a singular behavior of the 
SD band can be noticed also in the previous cranked HF calculation with the 
BKN force |TB|. In Fig. 6 we also plot the I =24 and 26 states, which are 



missing in Figs. 1-3, by smoothly extrapolating the / — tUrot and (/32,7) — / 
curves for the SD and the HD-like configurations, respectively (see, Ref. M 
for the treatment of the band-crossing region). 

The microscopic origin of this singular behavior may be understood when 
we examine the single-particle energy diagram in the rotating frame (routhians) 
presented in Fig. 7. We see that the rotation-aligned level associated with the 
[440] i orbit comes down in energy with increasing uirot and crosses the Fermi 
level at Urot — 2.9 MeVh which corresponds to / ~ 24. Thus, the yrast states 
above / ~ 24 are characterized by the occupation of the [440] | level by a single 
proton and a neutron. According to the deformed harmonic-oscillator model, 
N=Z=18 is a magic number associated with the HD shell structure with axis 
ratio 3:1, in which the [440]^ level is occupied by two protons and two neutrons. 
In order to distinguish the yrast states with / > 26 from the SD states below 
/ ~ 24 and keeping in mind a connection to the HD configuration, we call 
them "the HD-like configuration" , although the magnitude of the quadrupole 
deformation (32 obtained in the SHF calculation is in fact comparable to that 
of the SD shape rather than the HD shape. 

Let us note that if we regard the SD configuration as to correspond to the j- 
j-coupling shell model Ap-12h configuration 7r[(/7/2)^(sc?)~^] ® ^[{f7/2Y{sd)~^] 
(relative to ^°Ca) in the spherical limit, the maximum angular momentum 
that can be generated by aligning the single-particle angular momenta toward 
the direction of the rotation axis is 24^, and thus "the SD band termination" 
might be expected at this angular momentum. Interestingly, our calculation 
indicates that a crossover to the HD-like configuration takes place just at this 
region of the yrast line. 



5 Comparison with the CHO model 

The behavior at the high-spin limit of the SD band obtained in the SHF cal- 
culation possesses some similarities with that expected from the CHO model. 
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This model has been frequently used P5|-p0| as a simplified model of rotating 
mean fields. With obvious notations, the single-particle Hamiltonian of this 
model is written as 



1 

h' = hujk{c\ck + ^roth, (7) 
fc=i ^ 

where 

c\ = ^JmuJk/2h{xk - ipk/muJk), (8) 

with (xi,X2,X3) indicating {x,y,z), etc. 

The orbital angular momentum operator li consists of two parts: 

h = t''='^ + t''='^ (9) 

with 

/r^°) = ^^i7=(^3C.-4c3), (10) 

ir^'' = ^^1^(44 -CSC). (11) 

For a given value of Urot or Ih =< J2t=iih)i >, one can determine the os- 
cillator frequencies (cji, uj2, uj^) such that the selfconsistency condition between 
the density and the potential, 

AAA 
< E(^?)i >= ^2 < E(^2)i >= ^3 < T.(4)i >, (12) 

i=l i=l i=l 



is fulfilled under a volume conservation condition |5^. Here, the brackets de- 
note expectation values with respect to Slater determinantal states composed 
of single-particle eigenmodes of h'. 

Let us denote the total number of quanta in each of the three directions 
{k = 1, 2, 3) at tUrot = as 

^k=<J2i^ick + -h>, (13) 
i=i ^ 

and let us continuously follow the configuration specified by the set of values 
(Si, S2, S3) which are defined at uj^oi. 7^ as the number of quanta associated 
with the normal modes of the CHO Hamiltonian h! . In terms of S^, the 
selfconsistency condition at Wrot = is written as 

U^iSi = LlI'{L2 = (14) 

If the AA^ = 2 part of the angular momentum operator li is neglected, it is 
well known that there exists a maximum angular momentum = S3 — S2 for 
a given configuration (Si, S2, S3), where the shape is oblate and the symmetry 



axis coincides with the rotation axis This shape evolution is caused by 
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the effect of the AA^ = part of the cranking term, which tends to align 
the angular momentum of individual particles toward the rotation axis of the 
system (rotation alignment effect due to the Coriolis force). In the case of 
the doubly closed shell configuration for the SD magic number N = Z = 16 
(including the spin-degeneracy factor 2), corresponding to the SD band in ^^S, 
(Si, S2, S3) = (24, 24, 48) taking into account protons and neutrons. We would 
thus expect the "SD band termination" at the maximum angular momentum 
/c = S3 — S2 = 24. This number coincides with that evaluated in the previous 
section in relation to the j-j coupling shell-model configurations. 

On the other hand, the AA^ = 2 part stretches the system toward larger 
deformations, and actual shape evolutions as functions of angular momentum 
are determined by the competition and balance between these two effects. Fully 
taking into account both effects of the cranking term, Troudet and Arvieu 
5D[] found that there is a critical value ac of S3/S2, 



such that the configuration (Si, S2, S3) does not (does) reach the oblate limit 
if S3/S2 is greater (less) than Oc- This is because, for large deformations, 
the stretching effect of the AN = 2 term dominates at high spin over the 
alignment effect of the AN = term. In the case of ^^S, the SD configuration 
have S3/S2 = 002/^^3 = 2 > ftc at J = 0. Therefore, the "oblate limit" 
mentioned above will not be reached and the shape at the "band termination" 
point will be triaxial. 

Figure 8 shows the shape evolution in the {^2,1) plane, calculated for the 
SD configuration of ^^S in the CHO model. Here, the result of calculation 
for the configuration (22,24,54) is also presented, as an example of the HD- 
like configurations. We see that, although the triaxiality slowly sets in with 
increasing angular momentum, the shape of the SD states remains rather far 
from the oblate limit and exhibits a striking "back-bending" at about Ic = 24 
toward larger prolate deformations for / > Ic- Apparently, the behavior near 
the critical angular momentum 1^ for the SD band is quite similar to that of 
the SHF solutions presented in the previous section. On the other hand, it 
should be recalled in comparing Fig. 8 with Fig. 6 that the highest spin region 
of / = 26 ~ 32 on the yrast line corresponds to the HD-like configuration in 
the SHF solution: While the continuation of the SD configuration (24,24,48) 
to the / > 24 region as well as that of the HD-like configuration (22,24,54) to 
the / < 26 region are presented for the CHO model, only the yrast states were 
obtained and plotted in the SHF calculation. 

Figure 9 shows the angular momentum and the dynamical moment of in- 
ertia JT'^^-' as functions of the rotational frequency. We see that JT'*-^-' gradually 
decreases until the critical point. It is interesting to compare this property 
with that of J'^'^^ for the SD band in the SHF calculation (Fig. 5). Apparently, 
they are quite similar. This suggests that the gradual decrease with increas- 
ing cjrot of the dynamical moment of inertia for the SD band is rooted in the 
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existence of the critical angular momentum Ic associated with the quantum 
SD shell structure. We feel that a more detailed investigation of the SD states 
near the "band termination" point is a very important and challenging subject 
for a deeper understanding of the rotational motion of the nucleus as a finite 
Fermion system, 



6 Effects of time-odd components 

In this section we shortly discuss about the rotation-induced, time-odd compo- 
nents in the mean field. The moment of inertia of the SD band is expected to 
be a good physical quantity for identifying the effects of the time-odd compo- 
nents, since the pairing correlation plays only a minor role there. Concerning 
the effect of various time-odd components on the moment of inertia, we re- 



fer some recent works; for a semiclassical description, for a rotating 
nuclear matter, and for SD bands around ^^^Dy. 

Table 1 shows individual contributions from various kinds of time-odd term. 
It is interesting to note that the contributions from terms containing the spin- 
density, p{r), nearly cancel out with each other and, accordingly, the contri- 
bution from the current-density terms, denoted by B3 + B4, dominates in the 
sum. Such a remarkable cancellation of the spin-density terms was not seen in 
the case of ^^^Dy and may be specific to ^^S under consideration. 



In Fig. 10 we compare the results of calculation with and without the time- 
odd components. It is seen that the time-odd components increase the angular 
momentum for a given value of cUrot- Accordingly, the dynamical moment of 
inertia JT"*-^^ = rfJ/rfcUrot also increases. This trend is understood from the 
consideration of the local Galilean invariance of the Skyrme force |]5ll,^ (for 



a more general analysis not restricted to the Skyrme force, see p5| , |5^ ): If the 
time-odd components is neglected, the local Galilean invariance is violated and 
we obtain the moment of inertia associated with the effective mass m*. By 
including the time-odd components, however, the local Galilean invariance is 
restored and we get the moment of inertia associated with the nucleon mass m. 
The calculated result presented in Fig. 10 is consistent with this expectation, 
but a more quantitative analysis is not necessarily easy, because, as seen in 
Fig. 5, the calculated moment of inertia significantly deviates from the rigid- 
body value due to the shell effect. 
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Table 1 

Contributions from various terms in the time-odd energy density, 



B 



ll{Pv 



pD 



5i2pV + 5i3P"(p^ + P;), 



to the energy (in MeV), calculated at cj^ot = 1-0 MeV/h for the SIII and 
SkM* forces. Here, {jn,jp) and {p,^, pp) denote the nucleon currents and the 
spin densities (for neutrons and protons), respectively, and j = j„ + and 
p = Pn + Pp (see Ref. |Q for their explicit expressions). In columns designated 



by coefficients Bi, values after the spatial integration are listed, while the total 
value / drHoddi'f') and the sum of contributions from the current terms (the 
first two terms in the r.h.s. of the above equation) are shown in the columns 
denoted by "total" and "i?3 -|- B4" , respectively. For reference, the effective 
mass m* in nuclear matter for each force is also listed. 





B3 Bi Bq Bio Bu B12 Bn 


total 


B3 + Bi 


m* /m 


SIII 

SkM* 


-1.94 0.79 -0.17 -0.77 0.86 0.37 -0.18 
-1.83 0.90 -0.38 -0.44 2.45 0.00 -1.65 


-1.04 
-0.95 


-1.15 
-0.93 


0.76 
0.79 



7 ¥31 deformation 

As mentioned in Section 3, we found that the yrast states in the region 5 < 
/ < 13 possess a significant amount of non-axial octupole deformations of the 
lai type. It should be emphasized that such an exotic deformation is absent 
at / = but emerges at high spin. It has become possible to get this kind 
of solution by using the new cranked SHF code allowing for the simultaneous 
breaking of both axial and refiection symmetries. 



As in [39|, we define the octupole deformation parameters as 



A 

asm = iAn/3AR^) <J2ir'X3m)i>, (m = -3,---,3) (16) 

i=l 

with R = 1.2^1/3 Hg^g is a real basis of the spherical harmonics. 



^30 — ^30, ^3\m\ — "^C^S-lml + '^3-|m|)) 



1 

X3-\m\ = -^(X3\m\-Y^^^^), (17) 

where the quantization axis is chosen as the largest and smallest principal 
inertia axes for prolate and oblate solutions, respectively. The yrast solutions 
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in the region 5 < / < 13 have a^i 7^ but a^m = for m 7^ 1. (See |]5^ for a 
general discussion on this kind of deformation and its consequence on rotational 
spectra.) Figure 11 shows the calculated values of the Y31 deformation as a 
function of Umt- We see that the 031 value quickly rises when Urot exceeds 1 
MeV/h. 

The microscopic origin of the growth of the non-axial octupole deformation 
0:31 may be understood when we examine the single-particle energy diagram 
in the rotating frame (routhians) presented in Fig. 12. We note that a strong 
coupling and a quasi-level crossing between the rotation-aligned [330] | orbit 
and the [211]^ orbit take place near the Fermi surface in the region 1.0 < 
i^rot < 2.2 MeV/h. The matrix element of the r^Y^i operator between the 
two single-particle states is large, since they satisfy the selection rule for the 
asymptotic quantum numbers (AA = l,Anz = 2). This strong coupling is 
responsible for the a^i deformation appearing in this region of the yrast line. 

Figure 13 shows the potential energy function with respect to the 031 di- 
rection, calculated by means of the constrained HF procedure. Note the scale 
of the ordinate. Although we obtain a clear minimum at a finite value of 
Q!3i, the potential is rather shallow in this direction, so that the amplitude of 
the quantum-mechanical zero-point vibrational motion might be larger than 
the equilibrium deformation. If this is the treatment of dynamics go- 

ing beyond the mean-field approximation is required in order to investigate 
the consequence of the a^i deformation on the quantum spectra in the yrast 
region under consideration. This is beyond the scope of the present paper. 

It may be desirable to extend the potential energy curve in Fig. 13 to 
the ttai = limit. It turned out, however, difficult to do so, because many 
level-crossings take place with decreasing 0:31. [If we extrapolate to this limit 
assuming parabolic dependence on a^i, we obtain about 2 MeV as a very crude 
estimate of the energy gain due to the a^i deformation.] By the same reason, 
it is also difficult to follow the continuation of the Y^i band to the higher spin 
region as soon as it departs from the yrast line. 



8 Some remarks on experimental data 

Although rich experimental data are available for excited states of ^^S, the 
high-spin yrast region in which we are interested is rather poorly known at the 
present time. Accordingly, we discuss experimental references only briefly. 

For low-spin states with / < 7, detailed spectroscopic data are available up 
to excitation energy 11.76 MeV [|19],|2^. These excited states are shown to be 



well described by the spherical shell model calculations [19,0. In these works, 
some negative-parity states were interpreted as octupole-quadrupole phonon 
multiplets. As a matter of fact, we need to go beyond the simple mean field 
theory in order to discuss such spectroscopic data in the low-spin region. 

Highly excited states have been studied by various nuclear reactions as 
well as ^^O-^^O scatterings. Investigating the ^^0(2°Ne, a)^^S(a)2^Si (g.s.) 
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reaction, Morita et al suggested possible band structures of the quasi- 
molecular configuration of ^^O+^^O and of some parity-doublet-like structures 
with angular momenta 5~,6^, (7~), (8"*") at the 12-15 MeV region. Recently, 
Curtis et al ||56| investigated the region with 1=10-16 and the excitation energy 
32-38 MeV by means of the ^2c(24Mg,i60^^0)'^He reaction, and suggested an 
existence of highly deformed states in this region. It is tempting to compare 
these experimental data with our theoretical calculations. The experimentally 
explored regions are, however, about 10 MeV above the theoretical yrast line. 
Therefore, a more detailed spectroscopic study is needed in order to associate 
these data with the yrast structure. 



9 Conclusions 

We have investigated the high-spin yrast structure of ^^S by means of the 
cranked SHF method in the 3D Cartesian-mesh representation without impos- 
ing restrictions on spatial symmetries, and suggested that 

1) a crossover from the SD to the HD-like configurations takes place on the 
yrast line at angular momentum J ~ 24, which corresponds to the "band ter- 
mination" point in the CHO model, and 

2) non-axial octupole deformations of the V'si type play an important role in 
the yrast states in the range 5 < / < 13. 

In conclusion, we would like to stress again that the calculated yrast line 
for /= 14-20 lies about 10 MeV below the observed molecular resonance region 
associated with the ^^O-^^O configurations. Thus, a yrast 7-spectroscopy with 
higher resolving power is strongly desired in order to explore the high-spin 
region of the yrast line of ^^S . 
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Figure 1: (a) Excitation energy vs. angular- momentum plot for the yrast 
structure of ^^S, calculated with the SIII force. Density distributions on the 
plane perpendicular to the rotation axis are shown, as insets, for the SD band 
(solid line) and the Y31 band (dashed line). The calculation was done in step 
of Auj-ot = 0.2 MeV/h, and the calculated points (indicated by symbols) are 
smoothly interpolated by lines, (b) Same as (a), but with the SkM* force. 




Figure 2: Angular momentum I plotted as a function of rotational frequency 
cUj-ot for the SD band and the Y31 band in ■'^S. Results calculated with the SIII 
and SkM* forces are shown by solid and dashed lines, respectively. 
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Figure 3: Quadrupole deformation /32 plotted as a function of rotational 
frequency Urot for the SD band and the Y31 band in ^^S. Results calculated with 
the SIII and SkM* forces are shown by solid and dashed lines, respectively. 




Quadrupole Deformation 1^2 



Figure 4: Potential energy function at / = for the SD configuration (solid 
line) relative to that for the ground state configuration (dashed line) in ^^S, 
calculated with the SIII force. 
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Figure 5: Dynamical moment of inertia J''^'^^ = dl /duj^ot plotted as a function 
of LUrot for the SD band in ^^S. Results calculated with the SIII and SkM* 
forces are shown by sohd and dashed lines, respectively. For reference, the 
rigid moments of inertia Jrig — rn J p{r){y'^ -\- z^)dr with the calculated density 
p{r) are also indicated. 




Quadrupole Deformation ^2 



Figure 6: Shape evolution as a function of angular momentum, plotted in the 
(/92,7) plane for the SD and HD-like configurations in ^^S. Results calculated 
with the SIII and SkM* forces are separately shown. 
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Rotational Frequency Uiot [MeV/fi] 

Figure 7: Single-particle energy diagram (for neutrons) in the rotating frame 
for the SD band in "^^S, plotted as a function of cj^ot- The SIII force is used. Note 
that the structure of the yrast configuration drastically changes at Uj-ot — 2.9 
MeV/h, so that the diagram is discontinuous about this point, although levels 
characterized by the same asymptotic quantum numbers are linked by lines. 
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Figure 8: Shape evolutions as functions of angular momentum in the {(32, 7) 
plane, plotted with filled and open symbols, respectively, for the SD configu- 
ration (Ei,E2,S3) = (24,24,48) and the HD-like configuration (22,24,54) in 
the CHO model. 
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Figure 9: (a) Plot of angular momentum vs. rotational frequency in the CHO 
model. Solid line is used for the SD configuration (Si,S2,S3) = (24,24,48), 
while dashed line for the HD-like configuration (22,24,54). (b) Same as (a), 
but for dynamical moment of inertia JT'-^-' = dl/dcUrot- For reference, rigid 
moments of inertia — m < J2t=i{y'^ + > fo^ these configurations are 
also indicated. 
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Figure 10: (a) Angular momentum / plotted as a function of oUj-ot for the 
SD band in '^^S. Solid line with filled squares (dashed line with open squares) 
indicates the result with(without) the time-odd components. The SIII force is 
used, (b) Same as (a), but with the SkM* force. 
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Figure 11: Non- axial octupole deformation 0:31 plotted as a function of cjrot 
for the Y31 band in ^^S. Results calculated with the SIII and SkM* forces are 
shown by solid and dashed lines, respectively. 
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Figure 12: Single-particle energy diagram (for neutrons) in the rotating frame 
for the Ysi band in ^^S, plotted as a function of ujrot- The SIII force is used. 
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Figure 13: Potential energy function for the I31 band in ^^S at Urot — 2.0 
MeV/h, calculated by means of the constrained HF procedure with the SIII 
force. Note the scale of the ordinate. 
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